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The shear-modulus and yield-stress of amorphous solids are important material parameters, with 
the former determining the rate of increase of stress under external strain and the latter being the 
stress value at which the material flows in a plastic manner. It is therefore important to understand 
how these parameters can be related to the inter-particle potential. Here a scaling theory is presented 
such that given the inter-particle potential, the dependence of the yield stress and the shear modulus 
on the density of the solid can be predicted in the athermal limit. It is explained when such prediction 
is possible at all densities and when it is only applicable at high densities. These results open up 
exciting possibilities for designing in principle new materials with desirable mechanical properties. 



Introduction: Typical solids, whether amorphous 
or crystalline, respond linearly and elastically to small 
strains. Then the stress a%j increases like fxeij where \x is 
the shear modulus, is the strain 



conditions, to predict how fi and ay depend, say, on the 
density of the material. In addition we clarify the issue 
of universality of the ratio h/<jy which had come up in 
recent experiments on a family of metallic glasses [2]. 
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and u is the displacement field. Upon increasing the 
strain further the material yields and begins to flow; the 
average stress- value in the plastic steady state is denned 
here as the "yield-stress" cry. In crystalline materials 
this plastic flow is understood in terms of the motion 
of defects like dislocations [1]. It is much less clear in 
amorphous solids what the carriers of plasticity are and 
how to relate the onset of plasticity to the microscopic 
properties of the material. The aim of this Letter is to 
present a scaling theory that is able, subject to stated 



The yields stress and the energy drops. To ini- 
tiate the discussion we remind the readers of a typical 
strain-stress curve in amorphous materials. For this pur- 
pose we construct a glassy system consisting of poly- 
dispersed soft disks. We work with N point particles 
of equal mass m in two dimensions with pair-wise inter- 
action potentials. Each particle i is assigned an interac- 
tion parameter A^ from a normal distribution with mean 
(A) = 1. The variance is governed by the poly-dispersity 
parameter A = 15% where A 2 = ^tto^ ■ With the 
definition A y = + \j) the potential is constructed 
such as to minimize computation time [3] 
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where r y - is the separation between particles i and j. 
The shape of this potential for k = 8 and k — 10 is 
shown in Fig. 1 . Below the units of length, energy, mass 
and temperature are (er), e, m and e/ks where ks is 
Boltzmann's constant. The unit of time To is accordingly 
t = \J (m(cr) 2 /e). In the present simulations we choose 



Bq = 0.2 and employ an athermal quasi-static scheme 
[3] with Lees-Edwards boundary conditions. In Fig. 2 
we show a typical stress-strain curve for this model with 
N = 4096 and k = 10. For small strains the stress grows 
linearly with a slope this continues as long as the ma- 
terial is elastic, i.e. for stress less than ay. For stress 
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FIG. 1: Color online: The three different potential discussed 
in this Letter. 




FIG. 2: Typical stress-strain curves for a system with N 
4096 and k = 10 in Eq. (2). 



exceeds ay, plastic events begin to take place. In this 
regime we observe drops in the stress of size A CT , and we 
denote the mean size of these stress drops by (A a ). Be- 
tween the stress drops the strain is increased by intervals 
A e with a mean (A e ). Associated with the stress drops 
are potential energy drops which we denote by Au whose 
mean is (Au). In [3] the dependence of these averages 
on the system size was determined, with the results 



(A CT ) - CN? , (Au) = C v N a , 



(3) 



with C CT having the units of stress and Cjj of energy. 
The exponents satisfy a scaling relation a — (3 = 1 which 
can be easily derived from the exact energy conservation 
relation [3, 4] 



ay x (A e ) x V = ay x x V 
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Substituting Eq. (3) into Eq. (4), the experimental fact 
that ay and p are intensive quantities independent of N 
implies a — f3 = 1, and we can rewrite the last equation 
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FIG. 3: Color online: stress-strain curves averaged over 20 
independent runs for an athermal system with N = 4096, 
k = 8 (left panel) and k = 10 (right panel) as a function of 
the density, with the density increasing from bottom to top. 



in favor of the yield stress ay , 
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where p = N/V is the number density, Note that the 
combination fl^ 1 = pC\jjC a is a pure number whose 
non-universality will be discussed below. The challenge 
at this point is to provide a theory for the dependence 
of (i and the pre-factors C„ and Cjj on the experimental 
control parameters. In this Letter we focus on the ather- 
mal limit T — - and study the density dependence; the 
thermal theory is more involved and will be dealt with in 
a future publication. 

Results and Scaling Theory: In Fig. 3 we present 
the stress-strain curves of two systems, each with N — 
4096, as a function of the density, one with fc = 8 in 
Eq. (2) and one with fc = 10. We note the very strong 
dependence of both fi and ay on the density (60% change 
in the density results in a factor of up to 20 in ay). 
We also note the sensitivity to the repulsive law of the 
potential, changing fc by 25% changes the yield stress by 
a factor of 2. To understand this density dependence we 
point out that both ay and fi have the dimensions of 
stress, and from the potential the quantity having the 



same dimension is the function r 
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Thus if this 



function has a scaling behavior with r then ay and /i 
would have a related scaling behavior with p. If U(r) 
were a simple power law r~~ k this function would behave 



like r 
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potential Eq. (2) is not a simple power law but as we shall 
see the function r -1 9U q^ exhibits an effective power law 



of r G [p. 
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], see Fig. 4. We find that to a very 



or scaling behavior. We compute r 1 " U g r ' ' in the range 



FIG. 4: Color online: 
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[pmax , P m i„ ] ■ The line through the points represents the 



good approximation 
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with v = 4.80 for k = 8 and v = 5.87 for k = 10. With 
this effective scaling behavior, we can predict a-priori 
that being a number, should be independent of p 

and both cry and (/z) should scale like 



try 



(7) 



where (/i) denotes an average of the shear modulus over 
the elastic steps of the steady state plastic flow. Con- 
sequently we expect that the stress-strain curves would 
collapse by dividing all the stresses by p v . Fig. 5 demon- 
strates the perfect data collapse under this rescaling, to- 
gether with the high precision of the scaling laws (7). 

When scaling fails: To delineate the applicability of 
the scaling theory, we turn now to a third potential in 
which to the repulsive branch we add an attractive one, 
see Fig. 1. The potential chosen reads (with the same 
polydispersity of A y - ) 



U(r) = 



with U(rij) - 
2 1 / 6 A„- and 
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The attractive part U(r) 



We choose 
x 1 and the co- 



is glued smoothly to the repulsive part. 
U{r) = \P (^) where P(x) = Eto 4 
efficients Ai are chosen such that the potential is smooth 
up to second derivative [5]. We repeated the measure- 
ments of the stress-strain curves for this system, now with 
N = 2500, and display the results in Fig. 6, left panel. 
In the right upper panel we show what happens when we 
try to collapse the data by rescaling the stress by cry . Of 




FIG. 5: Color online: The same stress-strain curves as in Fig. 
3 but with the stress rescaled by p" , with v — 4.80 for k — 8 
and v = 5.87 for k = 10. The insets demonstrate the density 
dependence of cry and /i according to p v . 



course the stress-stain curves now all asymptote to the 
same value, but the curves fail to collapse, since \i does 
not in general scale in the same way as cry. The reason 
is the failing of the scaling hypothesis in this case: the 
calculation of r^ 1 ^^ in the range of r G [p m i x 2 ,p~[/ 2 ] 
fails to provide an effective power law, destroying the 
scaling hypothesis and the data collapse. Nevertheless 
even in the present case we can have predictive power for 
high densities. When the density increases the repulsive 
part of the potential (8) becomes increasingly more rel- 
evant, and the inner power law r~ 12 becomes dominant. 
We therefore expect that for higher densities scaling will 
be regained, and both cry and (/x) would depend on the 
density as p 7 . The two lower right panels in Fig. 6 show 
how well this prediction is realized also in the present 
case. 

We therefore conclude that given a potential one 
should be able to decide whether a scaling theory should 
be applicable and at what densities. Whenever the most 
probable inter-particle distance agrees with the average 
inter-particle distance scaling can be employed with im- 
punity. At low densities this condition may not apply 
due to the attractive part of the potential. Nevertheless 
with higher densities scaling and predictability of the me- 
chanical properties should get better and better. 

Finally we discuss the numerical value of the param- 
eter f2 = C a /(pCu)- We note that the data collapse 
indicates that this parameter is, to a very good approxi- 
mation, independent of p for a given potential function. 
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FIG. 6: Color online. Left panel: stress-strain curves for the 
potential (8) which has a repulsive and an attractive part. 
Right upper panel: demonstration of the failure of rescaling 
of the stress-strain curves. Lower panels: ay and {fj,} as a 
function of the density. Note that predictability is regained 
only for higher densities, the straight line is p 7 . 
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FIG. 7: Color online: The pure number f2 as a function of the 
density for the three potentials discussed in the text. Note 
that Q appears to increase with the exponent of the repulsive 
part of the potential whenever scaling prevails. 



For the two different potentials (2) with k = 8 and k = 10 
respectively, we find from our numerics that this parame- 
ter differs by about 5%, indicating non- universality. The 
lack of universality is even clearer with the last potential 
(8). In Fig. 7 we display the computed values of f2 for 



the three cases studied above. When scaling prevails the 
value of Q is constant up to numerical fluctuations. In 
the third case, where scaling fails, f2 is a strong function 
of p except at higher densities where scaling behavior 
is recaptured as explained. We can therefore conclude 
that the approximate constancy of Q found in a family 
of metallic glasses [2] , is not fundamental but only an in- 
dication of the similarity of the potentials for this family. 
In general f2 can depend on the inter-particle potential. 
It is quite clear from considering Eqs. (6) and (7) that 
the coefficients in the scaling laws (7) may well depend 
on the exponent k in the repulsive part of the poten- 
tial. The ratio of these pre-factors, being a pure number, 
could be independent of k, and f2 could be universal. It 
appears however that \x is increasing more with k than 
cry, and Q shows a clear increase upon increasing k. How 
to explain this dependence is beyond our present state of 
understanding and must remain an interesting riddle for 
future research. 

In summary, we have begun in this Letter to explore 
the relations between inter-particle potentials and the 
mechanical properties of amorphous solids that are made 
of these particles. Presently we focused on the athermal 
limit and showed how the dependence of the yields stress 
and the shear modulus on the density can be determined 
from first principles as long as the scaling hypothesis pre- 
vails. We also examined when this hypothesis failed, and 
pointed out that in any material the theoretically deter- 
mined relations can be trusted at high densities. In later 
publications we will examine the effect of temperature on 
the present issues. 
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